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Abstract 

The two-loop three-gluon vertex is calculated in an arbitrary covariant gauge, in 
the limit when two of the gluons are on the mass shell. The corresponding two- loop 
results for the ghost-gluon vertex are also obtained. It is shown that the results are 
consistent with the Ward-Slavnov-Taylor identities. 
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1 Introduction 



Calculation of radiative corrections to certain jet processes is becoming extremely im- 
portant (see, for instance, the reviews [HE] and references therein). In particular, the 
two-loop three-gluon vertex with two external gluons on shell (i.e. when two external mo- 
menta squared vanish) is a part of the next-to-next-to-leading order (NNLO) contributions 
to gg — > gg and qq gg processes (see, for instance. Fig. 1). The next-to-leading order 
(NLO) contributions to these processes were considered in ref. jSjl (see also in refs. ^). 
For a complete calculation of all relevant NNLO contributions, one also needs to calculate 
the two-loop four-point functions (see, for instance, in [SUED- 

Studying the three-gluon vertex (and other QCD vertices) in this on-shell limit is 
also important, in order to understand the structure of infrared (on-shell) singularities, 
to illustrate how the corresponding Ward-Slavnov-Taylor (WST) identities work in this 
divergent case and to develop calculational techniques that can then be generalized to 
deal with more complicated (e.g. off-shell) configurations. Performing the calculations 
in an arbitrary covariant gauge is useful, in order to exploit the consequences of gauge 
invariance (the WST identities) at all stages. In the calculation of physical quantities, 
the independence of the gauge parameter usually serves as an important check. 

The one-loop QCD vertices have been known for quite some time. The one-loop result 
for the three-gluon vertex, for off-shell gluons (with p\ = P2 = P3) and in an arbitrary 
covariant gauge, was presented in [71IH]. The general off-shell case, but restricted to the 
Feynman gauge, was considered in [Uj. Various on-shell results have also been given: in 
pUj . restricted to the infrared-singular parts only (in an arbitrary covariant gauge), and 
in ^1], with the finite parts for the case of two gluons being on-shell (in the Feynman 
gauge). The most general results, valid for arbitrary values of the space-time dimension 
and the covariant-gauge parameter, have been presented in an earlier paper ^21, where 
we have also collected the results for all on-shell limits of interest. Some results for the 
one- loop quark-gluon vertex (or its Abelian part, which is related to the QED vertex) 
can be found in [T8j . 

At the two-loop level, the QCD vertices were mainly studied in the zero-momentum 
limit ITK] . i.e. when one of the external momenta is zero. This limit is useful for 
studying the renormalization properties of QCD, since (for the considered vertices) it 
does not bring in any infrared (on-shell) singularities. In ^3], the renormalized results 
in the Feynman gauge were presented. In the unrenormalized and renormalized 
results for the three-gluon and ghost-gluon vertices have been obtained in an arbitrary 
covariant gauge, and the corresponding differential WST identity in QCD was analysed 
in detail. The relevant techniques for on-shell two-loop calculations have been studied in 
refs. |16| [T7| ITH] . In ref. ^H] , the two-loop electromagnetic quark form factor in massless 
QCD was considered. Corrected results for this form factor were later presented in [T^ I2Uj. 
In ref. [Ej , the two-loop scalar form factor was calculated. We also note that the structure 
and factorization properties of infrared singularities of two-loop order QCD amplitudes 
were recently discussed in ref. j21j . 

In the present paper, we discuss an algorithm to calculate two-loop three-point di- 
agrams with two external legs on shell. Then, we present the on-shell results for the 

^For a complete list of NLO contributions, see Figs. 6 and 8 of 
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three-gluon vertex and the ghost-gluon vertex, in an arbitrary covariant gauge, keeping 
the finite parts of the expansion in the dimensional regularization 22 parameter e. We 
consider the relevant WST identities in the on-shell limit and confirm that the results 
obtained are consistent with these identities. 

2 Preliminaries and WST identities 

The three-gluon vertex is defined as (see Fig. 2) 

rr4M3(pi'P2,P3) = -i 9 r^'^' r,,^,,,ip„p2,ps). (2.1) 

Here, the /"i^^aa q^j-q ^j-^g totally antisymmetric colour structures corresponding to the 
adjoint representation of the gauge group (for example, SU(A^) or any other semi-simple 
gauge group). Other colour structures do not appear in the perturbative calculation 
of QCD three-point vertices, at least at the one- and two-loop levels. To regulate the 
ultraviolet and infrared (on-shell) divergences occurring at the one- and two-loop levels, 
we shall use dimensional regularization 12^, with the space-time dimension n = 4 — 2e. 

Since the f°-^°'^°-'i are antisymmetric, ^2^3(^1 5^2,^3) must also be antisymmetric 
under any interchange of a pair of gluon momenta and the corresponding Lorentz indices. 
Therefore, in the limit of interest {pf = P2 = ^, P3 = P^) it can be presented as 

r/.iM2M3(Pl'P2,P3)|p2^p2^0,pi=p2 = ^MiM2 (j^l -P2)m3 ^l(P^) 

+ [^MiM3PlM2 - W3P2M1] t^2(p^) + [^miM3P2m2 - W3PIM1] Usip"^) 
+Plt,iP2^2iPl-P2)i,3 U^ip'^) + P2f,,Plf,,{pi -P2)t^3 U^{P^) 

+ [Pl^,^Pl^,2Pl^l■i - P2^Ji^P2^,2P2^l3] ^6(P^) + [PlmPlM2P2M3 - P2miP2a.2P1M3] U7{p^). (2.2) 

This decomposition is similar to eq. (29) of ref. |loj^ . All terms are explicitly antisym- 
metric with respect to (pi,/ii) ^ {p2,f^2)- At the lowest, "zero-loop" order, 

= 1, f/f ) = -2, [/f = -1, [/f = [/W = up = Uf^ = 0, (2.3) 

so we get the well-known tensor structure 

{P3-Pl) M2- 

(2.4) 

The lowest-order gluon propagator i^ 

^•"^(w.-«^). (2.5) 

where ^ = 1 — a is the gauge parameter corresponding to a general covariant gauge, 
defined in such a way that ^ = (a = 1) is the Feynman gauge. The gluon polarization 
operator is defined as 

n°r4(p)^- ^°^"^ {p'9,.,2 - p,.p,2) J{p'), (2.6) 

^For details, see section 4C and appendix F of ref. |12| . 

^Here and henceforth, a causal prescription is understood, ^/{p^ + iO)- 
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while the ghost self-energy is 

fiaia, ^p2^ ^ ^a,a, ^2 \^G{p^)] . (2.7) 

In the lowest-order approximation J^^^ = G^^^ = 1. 
The ghost-gluon vertex can be represented as 

f;:r"='(Pi,P2;P3) = -ig r"'"' r,^^{p,,p2;ps), (2.8) 

where pi is the out-ghost momentum, p2 is the in-ghost momentum, p^ and /is are the 
momentum and the Lorentz index of the gluon (all momenta are ingoing). For F^^g, we 
use following decomposition (see [q]!): 

At the "zero-loop" level, all the scalar functions involved in ()2.9j) vanish at this order, 
except one, a^^^ = 1. 

Whenever possible, we adopt the notation used in our previous papers jT3 E] • In 
particular, for a quantity X (e.g. any of the scalar functions contributing to the propa- 
gators or the vertices), we denote the zero- loop order contribution as (cf. eq. p.3|l ). 
the one-loop order contribution as X^^\ and the two-loop order contribution as X^^\ In 
this paper, as a rule, 

X(^) =X(^'«)+X(^'^), (2.10) 

where X^^'^^ denotes the contribution of gluon and ghost loops in a general covariant 
gauge ()2.5p (in particular, X^^'^^ corresponds to the Feynman gauge, ^ = 0), while X'^^''^^ 
represents the contribution of the quark loops. 

In general, the WST identity [SHI for the three-gluon vertex reads (see e.g. in ^^): 

r^i/.2M3(pi'P2,P3) = -J{pI) G{pI) [gf,|"''pl-Pl^,,pT] r^3M2(pi'P3;P2) 

+Apl)Gipl) [9,r'Pl-P2,.P'2] f^3m(P2,P3;Pi). (2.11) 

It is easy to see that the c and e functions from the ghost-gluon vertex ()2.9j) do not 
contribute to the WST identity 1)2.111) . In the on-shell case, some of the momenta squared 
vanish. Note that (in the case of massless quarks) J(0) = G'(O) = 1. In eq. 1)2.111) we can 
also consider permutations of the indices 1, 2, 3 corresponding to the contractions of the 
three-gluon vertex with or p2^ . To get all relations between the scalar functions in the 
case of interest {p\ = p\ = p\ = p"^), it is sufficient to consider the contractions with 
and Pi^ . 

The WST identity 1)2.11)) (contraction with pg^) yields just one condition on the scalar 
functions: 

U2{p') - Usip') + Ip'Ueip') + \p'U,{p') 

= -J(0) G{p') ^a(0,p^O) + V6(0,/,0) + Vd(0,p^O)l . (2.12) 



Adopting the notation used in ref . 0| , we have written in previous papers |12l I15| the arguments of 
these ghost-gluon scalar functions a, 6, c, d and e as momenta (vectors). In the present paper, it is more 
convenient to write the arguments as momenta squared. The order (3,2,1) of the arguments on the r.h.s. 
of eq. 1)2. 9|l corresponds to ref. ^ and has not been changed. 
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Considering contraction with , we get four relations: 

Uiip') + '^p'U,{p') = G(0) J{p') [a(0,0,p2) + lp%{0,0,p^) + '^p'd{0,0,p 

U2{p') = -2 G{0) J{p^) a(0,0,p2), 
U3{p')-h'Ur{p') = G{0)J{0) i/6(/,0,0)-G(0)J(/) \a{0,0,p^ 



(2.13) 
(2.14) 



i/rf(0,0,p2) 



P^Ueip^) = -G'(O) J(0) a(p2, 0, 0) - V c^(p^ 0, 0) 



+G(0) J(p') a(0,0,p') + ip2c?(o,0,p- 



(2.15) 



(2.16) 



Note that the l.h.s. of eq. ()2.12|) can be constructed from the Lh.s.'s of eqs. fl2.14|) - 
()2.16|) . This gives a condition on the scalar functions from the ghost-gluon vertex: 



G{p^) a(0,p^O) + y\0,p\0) + Vrf(0,p2,0) 



G(0) a{p\ 0, 0) + y%ip^ 0, 0) - ip'rf(p^ 0, 0) 



(2.17) 



Therefore, the conditions ()2.13p - ()2.17|) may be considered as a set of independent corol- 
laries of the WST identities (EZH). 

At the one-loop order, the diagrams contributing to the three-gluon vertex, two-point 
functions and the ghost-gluon vertex are shown in Figs. 1, 2 and 3 of ref. jT2j. The 
one-loop expressions for the Ui functions are presented in Appendix F of ref. |12j^ for 
arbitrary values of the space-time dimension and of the covariant-gauge parameter ^. In 
ref. ^2); they were obtained as a limiting case of the general off-shell expressions. Of 
course, they may also be obtained by direct calculation, using results for the on-shell one- 
loop integrals collected in Appendix A of the present paper. The corresponding one-loop 
expressions for the ghost-gluon scalar functions are, for the on-shell limits of interest, 
collected in Appendix B. They can be obtained from the expressions for general momenta 
presented in Appendix D of ref. jT^j, or by direct calculation. The one-loop expressions 
for the two-point functions can be found, for instance, in ref. [22] (see also in ^21 and 
Appendix C of the present paper). Collecting all the mentioned one-loop results, we have 
checked that they satisfy the conditions ()2.13|) - ()2.17|1 . in any space-time dimension n and 
for any ^. 

In the one-loop expressions, the following notation is use ci 



/t(p2) 



(n-3)(n-4) 



2N(n-4)/2 



e{l-2e) 



-p 



(2.18) 



-1) 



r(r2-3) 



r(3-f 



-■ye 



1 

12' 



^2^2 

vr e 



^' r(i 

7 



-2e) 
47 

1440' 



■vr e 



(2.19) 



^The divergent parts were presented earlier in ref. jlOj . whereas the Feynman-gauge resuhs, including 
the finite terms in e, are available in |1H . 

^Below, we shall also use the factor in two- loop results (they are proportional to 77^). Since the 
two-loop contributions involve poles up to l/e"*, we need the expansion of 77 up to the e'^ term. 
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where 7 ~ 0.57721566... is the Euler constant, whilst Cs = C(3) = ET=ij~^ - 1.2020569... 
is the value of Riemann's zeta function. 

We also use the standard notations Ca for the eigenvalue of the quadratic Casimir 
operator in the adjoint representation, 

jacdjbcd ^ ^ab _ ^ ^j^g gu(]v) group), (2.20) 

and Cp for the eigenvalue of the quadratic Casimir operator in the fundamental repre- 
sentation. For the SU(A^) group, Cp = {N'^ — 1)/{2N). Furthermore, 

T ^ NfTn, Tn = | Tr(/) = i , (2.21) 

where Nf is the number of quark flavours, and I is the "unity" in the space of Dirac 
matrices (we assume that Tr(J) = 4). 



3 Planar two-loop three-point integrals 

To calculate two- loop contributions to the three-gluon vertex (shown in Fig. 1 of ref. |15j). 
we consider contractions of eq. (j2.2|) with all possible tensor structures carrying three 
Lorentz indices /ii,/i2 and fi^. Note that non-planar graphs do not contribute to the 
two-loop vertex, since their over-all colour factors vanish, due to the Jacobi identity (see 
Fig. 6 of ref. |j22j, where this is explained). 

Technically, the problem is therefore reduced to the calculation of scalar integrals cor- 
responding to the planar two-loop vertex graph shown in Fig. 3a (and similar graphs with 
cyclic permutation of external momenta pi, p2 and ps). However, as a result of contract- 
ing the tensor structures, we get in the numerator some polynomials in scalar products 
of external and loop momenta. The complete basis for expanding these polynomials (see 
ref. |2Z|) includes (i) three external momentum invariants (e.g. pf, p\ and p|), (ii) six 
squared momenta corresponding to the six denominators shown in Fig. 3a, and (iii) one 
additional invariant, which can be chosen as q^. Diagrammatically, the latter member 
of the basis can be associated with the seventh line of an auxiliary "forward-scattering" 
four-point diagram shown in Fig. 3b. 

Since is missing in the original set of denominators (Fig. 3a), it always remains 
as a numerator, which cannot be cancelled against any of the denominators involved. 
Therefore, it is referred to as an irreducible numerator (see, for instance, ref. [2Z1)- In 
general, integrals with irreducible numerators require a special consideration 1^ . 
However, as we shall see below, this problem is not so serious in the on-shell case as in 
the general off-shell case, since the relevant "boundary" integrals can be calculated for 
any (integer) powers of the numerator q^. 

In terms of an algorithm, it is convenient to consider q^ as an extra denominator, 
remembering that its power z/7 is usually non-positive. Thus, let us consider integrals 
corresponding to the auxiliary diagram in Fig. 3b: 



K^iiji] z/i, z/2, 1/3, 1/4, 1^5, t'e, T^i) 



[{pi + rYY^[{pi + qYY^[{p2 - r)2]'^3[(p2 - qYY^{r-'^Y'^[{q - rYY^^q^y-^' 
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(3.1) 



where n = 4 — 2^ is the space-time dimension. We shall also need diagrams corresponding 
to the permutations of the external momenta in the diagram ()3.1|) . They are: 



K2{n;ui,U2,U3,U4,U5,U(i,UY) = Eq. (jSUl) with (pi,P2,P3) ^ (^3,^1,^2), (3.2) 
ifi(n; i/i, 1/2, 1^3, z^4, z^5, t'e, 1^7) = Eq. (jSH} with (^1,^2,^3) (P2,P3,Pi)- (3.3) 

The corresponding diagrams are the same as given in Fig. 3; the only thing to do is to 
permute the external momenta pi. Note that for the integrals K2 and Ki not all external 
lines of the four-point function are on shell, since some of them carry with p| = 7^ 0. 
Instead, one of the corresponding Mandelstam variables vanishes. 

To construct a procedure of calculating the integrals with different integer i/j, the 
integration-by-parts procedure [23 is usefulll If we introduce the notation j"*" to denote 
an increase of uj by one unit (and similarly let j~ denote a decrease of Uj by one unit), 
then the set of independent integration-bv-parts relations for the functions (jH.lj) can (for 
arbitrary momenta) be written as foUowqj: 



2/3! 3+ -Z^sl (2 -1 ) +(?7,-2z/i-Z/3-Z/5-Z/6 



P3Z/1I+ +pIu55' 



pluil^ +pIu33' 



-1/53-5^ 



-Z/6 



6+ (4- 



3 

5" 



+ Z/55+ 

1 - 2 



+ (n-z/i-2z/3-z/5-z/6) 
+ (n-z/i-z/3-2z/5-z/6) 
+ (n-z/i-z/3-f5-2z/6) 

Z/72^ 7^ + (n — 2U2 — U4^ — Uq~ Uj) 



P3I/22+ +P2i^77+ -z/22^4 +Z/66+ (^3 -4 ) -z/74 7+ +(ra-i/2-2z/4-z/6-z/7) 



i/22^ 



+ va4+ (3- 



Z/77^ 



+ (n-Z/2-Z/4-2z/6-Z/7) 



i^3 = 0, 

(3.4) 

i^3 = 0, 

(3.5) 

i^3 = 0, 

(3.6) 

i^3 = 0, 

(3.7) 

^3 = 0, 

(3.8) 

i^3 = 0, 

(3.9) 

i^3 = 0, 

(3.10) 

^3 = 0. 

(3.11) 

Analogous relations for K2 and Ki can be obtained by permuting the subscripts of pf^ 
according to eqs. ()3.2|1 and ()3.3|1 . 

Now, if we recall that we are dealing with the on-shell case, pi = p'2 = 0, we can see 
that some terms on the r.h.s. of eqs. ()3.4j) - ()3.1H) (and in the analogous relations for the 
K2 and Ki integrals) vanish. In this case, the following symmetry property is valid: 

Ki{n] ui, 1/2, z/3, z/4, z/5, z/6, z/7) = K2{n; z/3, z/4, z/i, 1/2, 1/5, i/e, z^?) (when = pj). (3.12) 

Therefore, it is sufficient to consider the and K2 integrals. 

^The application of this procedure to the calculation of the integral Ks{n; 1,1,1,1,1,1,0) was presented 
in JH| (see also eq. HH.13|I below). In a more general context, it was also discussed in [M)] . 

^The set of recurrence relations for calculating one-loop three-point functions was considered in 
ref. [SI- 



pIu22^ +pIva4.-^ -z/22+7- -1/44+7 +Z/66+ (5- - 7") +(n-z/2-z/4-z/6-2z/7) 
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When certain powers of propagators Vi vanislJl, the corresponding integrals Ki can 
be calculated in terms of F functions or finite sums over terms involving F functions. A 
collection of relevant results for such "boundary" integrals is presented in Appendix D. 
Using the relations ()H.4j) - ()H.lH) . the integrals Ki with integer powers of propagators z/, can 
be reduced to a set of such "boundary" integrals, or to analogous integrals where some 
z/'s are negative, i.e. the corresponding denominators are in the numerator. The latter 
integrals can also be reduced to boundary integrals (with the corresponding i/'s equal to 
zero), by using the tensor decomposition in appropriate self-energy-type sub- loops, i.e. 
formulae similar to those collected in Appendix A of ref . ^Hj (see also |29t l32] ) . 

It should be noted that, starting from integrals with non-positive z/7 and using relations 
()3.4|1 - ()H.11|1 . we never get integrals with positive z/7, since 7"*" is always accompanied by 
z/7 and one cannot "overcome" z/7 = 0. Nevertheless, in some cases integrals with positive 
z/7 may appear, due to the above-mentioned tensor decomposition in the sub-loops. Since 
the relevant boundary integrals can be calculated for an arbitrary z/7 (see Appendix D), 
this does not create extra problems, even if all seven z/'s are positive (see eq. ()D.6|l ). 

Let us illustrate this procedure by some important examples. One of them is the 
integral ^'3(1, 1, 1, 1, 1, 1, 0) (Fig. 3a), which was calculated in [121 IT7| ITH]: 



i^3(n; 1,1, 1,1, 1,1,0) 
1 

2e 



— [K^{n- 0, 2, 1, 0, 1, 2, 0) - 2K^{n- 0, 1, 0, 2, 2, 1, 0) + 2K^{n- 1, 1, 2, 2, 0, 0, 0)] 



. (3.13) 



To reduce the number of terms, some obvious symmetry properties of the boundary 
integrals have been used. Another example is the K2 integral with the same z/'s: 

ir2(ri; 1,1, 1,1, 1,1,0) 

^ [i^2(n; 1, 0, 0, 2, 1, 2, 0) + K2in; 0, 1, 0, 1, 2, 2, 0) + 2K2(n; 1, 0, 0, 3, 1, 1, 0)] 



2^(1 + 25) 



3 3 3 In ^ Iq ^ ^ / \ 



. (3.14) 



It corresponds to Fig. 3a, with the left lower momentum off shell. Expanding the t] 
factor ()2.19p . one can see that our results ()3.13p - ()3.14|) coincide with those presented in 
[T6| \T7\ HHj. Our eq. ()3.13|1 corresponds to the result for diagram 6A (with numerator 
= 1) presented in the Appendix of ref. [16 , to the first line of Table 4 of ref. ("fig- 2"), 
and to eq. (12) of ref. ^H]- Our eq. (j3.14j) corresponds to the second line of Table 4 of 
ref. [m ("fig. 4"), and to eq. (27) of ref. [180. We have also checked the results for 
other scalar integrals listed in refs. [T7\ [T8j , namely: all results for the diagrams from 
6A to 3A in the Appendix of ref. (including those with numerator^], but excluding 

^Diagrammatically, this can be understood as shrinking the corresponding line to a point. 

-'^'^In eq. (27) of ref. 18,, 9C4 should read 9C2- We also note another obvious misprint: in the last diagram 
on the r.h.s. of eq. (9) of ^H|i the "upper" line should contain a dot, i.e. the power of the propagator is 
equal to 2, which is clear from their eqs. (8) and (10). 

^^We note a misprint in the result for the diagram 4C with numerator (I ■ p), 15/(16tj) should read 
5/(16aj) (in our notation, lo ^ In addition, P is forgotten in the product of denominators in 

diagram 5F, and the denominator [I — r)'^ should read (Z + r)^ in diagram 4D. 
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the non-planar diagram 6B); the three remaining hnes in Table 4 of ref. j^, as well as 
the results listed in Appendix C of PTlpI (excluding the non-planar stuff, eqs. (C.6) and 
(C.7)); eqs. (26) and (28) of ref . [IB]. 

The algorithm discussed in this section makes it possible to evaluate all relevant inte- 
grals (corresponding to planar two-loop diagrams) for an arbitrary value of the space-time 
dimension n. However, since in physical applications one usually needs the expansion in 
e = (4 — n)/2, we present below the results for the two-loop vertices in an expanded form, 
up to i.e. keeping the poles and the finite terms. 



4 Two-loop results for the three-gluon vertex 

Below, we list the unrenormalized two-loop contributions to the functions UiijP') occur- 
ring in the three-gluon vertex. They were calculated using a set of REDUCE [231 P'^O" 
grams based on the algorithm described in the previous section. The two-loop diagrams 
contributing to the three-gluon vertex are shown in Fig. 1 of ref. ^Sl- The results are 
expanded in e up to the finite terms. The factor rj is defined (and its expansion in e is 
given) in eq. ()2.19|) . The colour factors Ca, T and Cp are defined at the end of section 2. 

The contributions of the diagrams without quark loops, in an arbitrary covariant gauge 
(for the definition of the gauge parameter ^, see eq. fj2.5|) ) are 

rr(2£), 2N ^2 9^ ^ 2n 2rfl /21 9 A 1 /II 335, 25,. 



(47r)"' ' V64 128 V Vl2 384^ 256 

1 / 37 19 2 623 3 ^ 103 2 1 
72 V 72 "^192^ 576^ ^ 128^ ^ ^ 128^ ^32^ 



1 / 239 1 2 11. 3343, 37 2. 

+- \ Tx H Cs H ^ H TT^f 

£ V 54 72 4 ^'^ 1728^ 384 ^ 

i^/- f _|_ ^37 2 + —f'^ 

32^^^ 384^ 256 ^ 16^ 

6653 119 2 139. 21 4 1333 . 37 191 



324 432 12 " 320 5184" 192 " 32 
64 ^ 1152^ 128 ^ 64^'^^ 16^ 16^ J ^ ^ ' 



- ^^(47r)"^ 1 166^ ^3 1,16 + 8^; 



1/15 5 2 1 . 1 2f 61 2 



(47r)" 

■{-■ 

V 8 

1 /593 1 o 7. 67 . 1 1 . . 3U^ 13 



£2 V 8 48 12^ 192 ^ 64^ 



£ V 48 36 4^^ 288^ 12 ^ 32^'^^ 48^ 16^ 

17939 19 2 107, 1 4 8381, 31 2. 

H \ Ti Cs H ^ ttV - 5C3£ 

288 108 12 20 1728^ 96 ^ 

^ -y^ - - ^ye + ^M' - ^ie - \e] + o{e), (4.2) 



640 " 288" 192 " 16 
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In eq. (C.5) of ref. 17 , in the second term in the square brackets (the term containing D) the 



numerator 2 + should read 2 + (note that the e used in ^7] corresponds to our — 2e) 
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1 /79 65 . 103 27 



£3 V96 192^ 192^ 256^ J 



1 + 11^2 , 157 13_ , _ 75]_ , _ J_ _ ^.3 _ 1.4 
.72 48 288^ 128 ^ 1152^ 192 ^ 256^ 32^ 
1 / 6125 155 ^ 13 . 97 . 1 135 



6 V 864 72 4 ' 1728" 48 " 64 



864 " 128 " 32'" 32" 768 

265205 206 o 227, 1 . 140815, 319 239,, 



5184 27 12 " 8 10368 288 " 32 

71 4. _ 23687 1 _ 1 . ,2 _ J_ 4,2 _ 145 3 

1280 ^ 5184 ^ 6 ^ 64^^^ 640 ^ 64 ^ 

5 o . q 5 , . q 13 



768 



„2..(2,^).2^ _ ^2 9^V\ „2w2. f 1 , , , 55 ^ 5 

^ ^ - (47r)"^ ^ ^ l£4 I32 + 128^ + 128^ ) 

1 /II , 301 539 65 ^ A 
£3 1^16 ^ 96 ^ 192^ 48^ 128^ 7 
1 / 9 11 2 389, 9 2. 1541,2 49 2,2 
+i^lT6 + 96^ -144^+16^^+1152^ +192^^ 
4267,3 1 2.3 3 ,4 1 ,A 
2304^ 32 ^ 16 16 / 

1 /2195 271 . 19, 34549, 7 15,, 7031,. 5 

+- TT^H Cs ^+-7rVH C3^H ^ TT^^ 

e V 96 144 8 ^ 1728 ^ 9 ^ 2 ^'^^ 1728^ 576 ^ 

221 ^2 20881 3 1 2,3 21 . 293^4 1 2,4 3 A 
+"6r^^^ ^U56^ "16''^ +64^=^^ + 256^ + 384^^ + 32^ J 
84683 2473 2 §9 1 4 1212661 685 ^ 

+ 576 - ^16"^ + 24^^ + 16^ - 10368 ^ + 216^ ^ 
1661,, 19 4, 39857,2 1327 2^2 109, ,2 

+^^^^ + 80^ ^ - + 3456^ ^ + ^^'^ 

139 4,2 87085 3 73 2.3 21 ,3 3 4,3 

+ 1280^^ +^^ + 768^^ "32^^^ + 256^^ 

+i^^+Tk-^^^-|^3e^+l^1+^(^)' ^'-'^ 



(47r)nv ^ ' I £4 V8 64^7 £3 V24 192'' 128 
1/10 1 2 653 7 2, 67 2 1 3 
+ £2 It - 16^ - 288^ ^ 192^ ^ - 128^ ^ 16^ 
1 /385_ n 21 _4303 5 5 .^.2 , 1 2.2_1.3\ 

£ V216 72 8 ^^ 864 ^ 192 ^ 4^^^ 64^ 128 ^ 16^7 
9557 41 2 113. 9 4 12889^, 9 2, 103., 9 4, 



1296 108 12 "" 160 1296 " 32 " 16 " " 320 
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457 2 7 o.o 9 



192^ 128 



r^'e + j^C3e-e] + 0{e), (4.5) 

2 



P^^ - u^.S P ) 132^4^ +,3[ 3 + 192^ + 128^ ) 



+ 



1/157 1 2 53 9 2> 155 2 13 2^2 29 
i2 1^144 ~ 48^ ~ 144^ ^ 64^ ^ ~ 128^ ^ 384^ ^ ^ 256^ ) 



1 /2365 1 2 1^ 3565^ 5 2^ 33^^ 
+ llW"36^ -8^-^- 1728^+64^^+32^^^ 

~192^ + 128^ ^ + 64^^^ ~ 64^ 
40429 _ ^^2 , 31 _ J_^4_ 150367 269^2. _ 3 
5184 432 12^'^ 160 10368 ^ 288 2^^^ 640 ^ 



(47r)"' ' V32 32^ 16 

1 /193 139^ 289 . 27 . 

H — \ £ £ f 

£3 V 48 192^ 384^ 128^ 

1 /125 83 2 119 _ I 2, _ 355 2 2.2 _ 133 3 

144 96 144^ 24 ^ 576^ 384 ^ 256^ 16 
1 / 29137 311 2 19^ 2027^ 43 2^ 49^^ 
+ £1-^ + 72"^ + y - 1728^ + 192^ ^ - 16^^^ 
8621 .0 43 7 . .0 39., 1 



L6^ ; 



^ + 384^^ +64^^^ "64^ " 384^^ " 4^ J 
1100617 6377 2 337, 53 . 50779, 19 2. 



1728 



5184 432 6 ^ 160 10368^ 32 

107,, 19 4, 15529 2 55 2,2 1 . .2 7 ,,2 



16'" 320 " 1296 " 96 " 32'" 1280 

32^ 768 ^ 16^-^ 32^ J ^ ^ ^ ' 

The unrenormalized two-loop contributions of the diagrams involving quark loops are 



(47r)"' ' I £3 V24 6^/ £2 13 
,1 + 1^2 _ 85 2 \ _ 449 43_ 2 , 19 . _ 197 41 2 
e U08 72 54^ 3^ / 324 108 3 162^ 18^ 



, 1^2 ,5 4 2X _ 2113 _ 23_ , _ 85 190 44 1 
U2 36 9^ 3^7 72 108 6 ^^ 27 ^ 9 ^ J 

(47 
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1/35 1 2 1> 11>2A 1 /4405 151 2 . 23 47 2 



^e^\ 36 6^ ^ 18^^ 18^ y ^ £ V 216 72 ^ 108^^27^ 
152755 1661 o 133, 1 4 235 



1296 216 12 '" 20 324 

4 ^2 



(47r)"' ' U3 V4 24^ 12^ 3 
1 /5125 53 o , 457^ 17 1111 . 1 25 



1 /5125 53 2 457 17 2,, 1111,2 1 2,2 , ^5 

153277 3307 o 113. 1 4 3346 31 2. 28,, 

+ Cs ? + — 

144 108 6 10 81 ^216 3 ^ 
^ 1919 .0 1 0.0 7 _o 131 



1/64 1 2 100 A 8101 37 2 1^ 767, 1,2! 
+ 11-27-36^ +^V+W-54^ +3^^ + ^^-3^ I 



4 ^2 



-4Ci.r^(V)-" + 0(£), (4.12) 



1 /2387 13 2 49, 1 2. 7,2 

+ - TT^ H f n^f - -r 

e \ 216 18 27^ 24 ^ 6^ 

^42515 223 o 35. 1519. Ho... 77 



1296 108 

-.4 ^2 



2 35^ 1519, 11 2, ^ , 77,2! 
^ +T^^ + l62-^-72^^ + ^^^-18n 



(47r) 

1/37 1 2 167, 11,2 

+^(-36-3^ -72"^ + ¥^ 
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1 /9419 125 2 . 61 1 2. 107,2\ 

+- - 2C3 — i + —T^ i + — i 

e \ 216 36 ^ 54^ 24 ^ 54 ^ / 

334055 352 2 29 , 1 . 263^ 5 2. 457 . 

+ Ca + i + —T^i + SCa^ + C 

1296 27 6 ^ 10 81 ^ 24 ^ 162^ 

+C,T ^S-P'r- (4 - \^') + 54 - A.' - 52C3} + 0(e).(4.14) 

We have also obtained the two-loop results for the ghost-gluon scalar functions in all 
on-shell limits of interest. They are given in Appendix E, whereas the relevant two-loop 
contributions to the two-point functions JijP') and GijP') are collected in Appendix C. 
Using all these expressions, together with the one-loop contributions, we have checked 
that all the results obtained satisfy the WST identities ()2.13p - ()2.17p . as they should. 

The renormalization of the results for the three-gluon vertex (and other three- and 
two-point functions involved) was discussed in detail in section 8 of ref. The cor- 

responding renormalization factors (Zi, Zi, Z3 and Z^) in the MS (or MS) scheme have 
been presented in refs. jHH ESj- For a detailed discussion of these results, together with 
a list of misprints, see Appendix B of ref. ^3] . To construct renormalized expressions for 
the three-gluon vertex functions f/j at the two-loop level, we need 

(i) to take the sum of (unrenormalized) zero-, one- and two-loop contribution^^ consid- 
ering the coupling constant and the gauge parameter as "bare" quantities, g ^ Qb and 

(ii) to substitute and in terms of the renormahzed g and ^, multiplied by the ap- 
propriate Z'-factors (see eqs. (8.8) and (8.9) of ref. [T3|): 

(iii) to multiply the resulting expression by the corresponding Z-factoi0, namely Zi (see 
eq. (B.l) of ref. [E]). 

Since the resulting renormalized expressions are as cumbersome as the unrenormalized 
ones (and can easily be obtained from the latter ones), we do not present them here. They 
also contain infrared (on-shell) poles in e up to l/e'^. 



5 Conclusion 

In the limit when two external gluons are on shell, we have calculated the two-loop 
contributions to the three-gluon vertex, in an arbitrary covariant gauge, keeping finite 
terms of the expansion in e = (4 — n)/2. In this limit, the three-gluon vertex is described 
by seven scalar functions f/j(p^) associated with different tensor structures, see eq. ()2.2|) . 
The results (listed in section 4) contain on-shell singularities up to The ultraviolet 

singularities are at most and should be removed by the renormalization. In a realistic 
physical calculation of squared amplitudes, the infrared (on-shell) singularities should be 
cancelled by the contributions of one-loop diagrams with soft emission from the external 
legs, etc., according to the Kinoshita-Lee-Nauenberg mechanism [30]. In this way, our 
result will be useful as a "block" in the calculation of NNLO corrections to physical 
amplitudes. 

^^Note that the one-loop expressions should be expanded up to terms, since they may be multiplied 
by other one-loop contributions involving poles. 

"'^^For the ghost-gluon vertex (see Appendix E), the renormalization factor Zi is required (see, for 
instance, eq. (B.2) in ref. |15p. 
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We have also calculated the ghost-gluon vertex ()2.9|) in all on-shell limits of interest; 
the results are collected in Appendix E. We have confirmed that the obtained results obey 
the corresponding WST identities ()2.13p - ()2.17p . 

We note that in the on-shell case considered, the problem of irreducible numerators 
in three-point two-loop integrals already shows up, but it can be overcome in a relatively 
simple way, since the relevant boundary integrals can be calculated for any integer powers 
of this numerator (see Appendix D). In the zero-momentum calculation p!3], there was no 
such problem at all. In the general off-shell calculation, though, the problem of irreducible 
numerators is much more severe |27j . 

Our results can be considered as a further step, in addition to [El CHI, towards cal- 
culating the two-loop QCD vertices in more complicated cases, such as the on-shell limit 
with just one gluon on shell, or the general off-shell casJ^ 
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-'^^In principle, the techniques for calculating the relevant off-shell scalar integrals are already available 

123 EZI 
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Appendix A: Some useful one-loop formulae 

The "triangle" integral with the external momenta pi, p2 and p^ = —pi — p2 is defined as 



J{n; z/i, z/2, z/3) 



d"g 



[iP2 - q)Y'[iPi + QyPiQ^)"'' 



(A.l) 



where n = 4 — 2^ is the space-time dimension. For general values of n and Ui, the result for 
the integral ()A.1|) can be expressed in terms of hypergeometric functions of two variables 

When two external legs are on shell, Pi=P2 = (Pa = p"^), the following simple 
formula can be easily obtained: 



J(ra;z/i,z/2,z/3) 



pI=pI=o 



jl-n^n/2(^p2^|n/2-S 



.:r(t-^l-^3)r 



-Z/2-Z/3 



r(z/i)r(i.2)r(n-Ej^^) 



(A.2) 

In particular, we shall need this formula for the case when one of the indices is a negative 
integer, z/3 = —s. This means that {q"^)^ is in the numerator. When 1/3 = 0, the r.h.s. of 
()A.2jl gives the well-known result for the one-loop two-point function. When z/i or i>2 are 
non-positive integers, we get zero. 

We also need the result for the triangle integral with one leg on shell. Assuming that 
pI 0, pI ^ and pi = 0, we get 



J(n; z/i, z/2, Us] 

r(t- 



Vl - 



Jl-r^^r^/2J^^2^Jn/2-EI., 



r t-z/i-z/2 



r(E 



^i-2 



T{V2 

r [vi + z/3 - 



z/i, Sz/j - I 

Z^l + - t + 1 

f)r(f-z/3) / 



Pi. 



J^2, 2 ~ 
H - Z/i - Z/3 + 1 



Pi, 



(A.3) 



where 2-^1 is the Gauss hypergeometric function. If we use the well-known formula of 
analytic continuation of 2-^1 function from the argument z to 1 — z (with z = P2/P1), 
we reproduce the result presented in Appendix A of ref. [Hn] (namely, their i^-integral at 
= 0). When z/i = —s (where s is a non-negative integer), the second term in the braces 
of ()A.3|) vanishes, and we get 



J{n; -s, z/2, z/3) 



r (f + s - Z/2) r (f + s - Z/3) r (z/2 + Z/3 - s - f 



r(z/2)r(z/3)r (n + s - z/2 - z/3) 

with a terminating 2-^1 series containing {s + 1) terms. 



-S, Z/2-hZ/3-S-: 
Z/3 - S - f + 1 
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Appendix B: One- loop results for the ghost-gluon 

vertex 



At the zero-loop level, we have 

flW (p2^ 0, 0) = a(°) (0, 0) = (0, 0, 



(B.l) 



whereas all other ghost-gluon functions are equal to zero at this order. 

The diagrams contributing to the ghost-gluon vertex at the one-loop level are shown in 
Fig. 3 of ref. ^21- The general expressions listed in Appendix D of ^2] give the following 
results in the on-shell limits of interest: 



a«(/,0,0) 



5(1) 
5(1) 
5(1) 



0,p2,0) 
0,0,/) 
p',0,0) 

0,p^0) 

0,0, p') 
p',0,0) 
0,p\0) 

0,0, p') 
p\0,0) 
0,p\0) 



d'^'\0,0,p') 



(4vr)"/2 16(n-4) 



X 



4(n - 4) - 2^{2n^ - On + 8) + ^'^{n -2){n- 4) 



^2 r] 



Ca 



(47r)"/2 4(n-4) 
Ca 



k{p^) 2 + ^(n^ - 6n + 7) 



(47r)"/2 8(n-4) 
Ca 



g"^ T] 



(47r)"/2 16(n-4)p2 



[2(3n-8)-e(n-4)] 



16(n - 4) + 4^(2^2 - 15n + 30) - ^^(n^ - 8n + 20) 



^2 



(47r)"/2 8(n-4)p^ 



g^ T] 



n-A)+ 4^(n^ - 7n + 13) - ^^(n^ - 8n + 14) 
Ca 



(47r)"/2 4(n - 4)p2 
Ca 



g^ T] 



(47r)"/2 4(n-4)p2 
Ca 



K^p") [s + 2i{2n-h)-e 
/t(/)e(n-6 + 20, 



^2 r] 



(47r)"/2 8(n-4)p2 



4(n - 6) + 2S,{n - 2)(n - 6) + ^^(n^ - lOn + 20) 
Ca 



(B.2) 
(B.3) 

(B.4) 

(B.5) 

(B.6) 
(B.7) 

(B.8) 
(B.9) 



^2 n 



(47r)"/2 16(n-4)p2 
Ca 



K{p^){n-Q) 8 + 4^(n-2) +^^(n-4) ,(B.10) 



g'^ T] 



(47r)"/2 16(n-4)p2 
Ca 



k{p^) 8(r2-6)-4^(5n-18) + r(^ -4n-4) (B.ll) 



g^ T] 



(47r)n/2 8(n-4)p2 



8(n - 6) + 2^(2n^ - 15n + 32) - C{n - 8n + 14) 
/s:(p2) [2(n -6) -2^(2/1-5) +r 



(4vr)"/2 4(r^ _ 4)p^ 



(B.12) 
(B.13) 
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e(^)(0,p^O) 



Ca 



(47r)"/2 2(n-4)p2 



nip") 2 + e(2n - 9) + 2^' 



P 



(B.14) 
(B.15) 



where i] and k(p^) are defined by eqs. ()2.19|1 and ()2.18j) . respectively. These results are 
valid for arbitrary values of n and ^. Note that there are no quark-loop contributions at 
the one-loop level. 



Appendix C: Two-point functions 

For arbitrary values of n and ^, the results for the one-loop two-point functions (see 
eqs. ()2.(ij) and ()2.7j) ) are available elsewhere (see |3Sl ESI IE])- For completeness, we also 
present them here: 



'''' <P'){-^ 



(47r)"/2 



4(3n-2) 



n 



+ 4(2n-7)^- (n-4)^2 



K{p') [2 + {n- 3)^] , 



+ 2T 



n — 1 



(C.l) 
(C.2) 



(4vr)"/2 4 

where nijp') and rj are defined in eqs. ()2.18|1 and ()2.19|1 . respectively. The corresponding 
diagrams are shown in Fig. 2 of ref. |12j . 

Two-loop diagrams contributing to the gluon polarization operator are shown in Fig. 3 
of ref. Calculating their sum, we get the following unrenormalized results [T3] : 



583 113^ 19 . 3 . 

\ f ^ 

72 144^ 24^ 8^ 



14311 
432 



1961 



1 /lOl 8 
~e VTs" ^ 9 
142 



2 

3' 

22 



27 



9 



(C.4) 



where C^i, T and are defined at the end of section 2. 

The two-loop ghost self-energy diagrams are shown in Fig. 4 of ref. ^Sl- Their sum 
yields the following unrenormalized results: 



4 ^2 



2 9_V_ 
(47r) 



1 /5 7 



1 



599 
"32" 



9 



4" 64 
17 



32 



1 /83 7 
e VI6 ^ 32" 



16 



rC3 + 0{e) 



(C.5) 



G'^-Hi-) = CUT^ (-i-)-^- _ ^ _ 1} + (C.6) 

The discussion of renormalization, as well as the renormalized results for J and G can 
be found in section 8 of ref . (see also [Hj , where these results in the Feynman gauge 
are presented). 

Appendix D: Boundary integrals 

Below the case Pi = P2 = 0, p| = is understood. Using the integration-by-parts 
relations ()3.4|) - ()3.11|) . we can reduce the integrals with six denominators to integrals 
where some of the lines are shrunk, i.e. some of the z/'s vanish. Here we list explicit 
results for such "boundary" integrals of interest. Apart from the one-loop formulae listed 
in Appendix A, the following well-known (see, for example, in |18| IH]) trick is useful, to 
combine two internal lines attached to a triple vertex with the remaining external line on 
shell: 



1 



1 



r(z/ + z/') } a^-Hl-aY'-^ da 



^) r('^') i [{q - apYr 



In fact, this is nothing but the Feynman parametrization of a product of two propagators, 
where it is taken into account that the difference of their momenta is light-like. 

The diagrams formally corresponding to the boundary cases of the integrals and 
K2 (listed below) are drawn in Fig. 4 and Fig. 5, respectively. To distinguish an off-shell 
external line (corresponding to ^3) from the on-shell ones (corresponding to pi and ^2), 
the former is drawn as a double line, which can be associated with the sum of pi and p2- 
Note that when using the tensor decomposition in the sub-loops we may get some integrals 
with positive v-j. This is not a problem, because the relevant integrals can be calculated 
for any v-j (see below). When the result is valid for an arbitrary z/7, the corresponding 
line is solid; when it is valid only for non-positive integer z/7, this is indicated by a dashed 
line, as in Fig. 3b. 

The results (ID.2j) . ()D.4jl and ()D.9jl can be obtained by repeated use of the one-loop 
formulae (see Appendix A). To get other results, the trick (jD.lll has been used. For the 
integrals ()D.3|) . ()D.7|) . ()D.8|) and ()D.10|) . this trick was used for two pairs of propagators. 
The notation Ei/j means I]J=i the sum over all i/'s involved (excluding those equal 

to zero). In particular, if z/7 = —s then Sz/j = Ylt=i — s. 

D.l: K-^ integrals 

The following boundary integrals can be expressed in terms of one-term products of F 
functions: 

K:,{n- z/i,z/2,z/3,z/4,z/5,0,z/7) =i'-'"vr"(p2)"-s->r(z/i + z/3 + z/5-t)r(z/2 + z^4 + ^^7-t 

r(f-z/i-z/5) r(f -Z/3-Z/5) r(f -z/2- z/7) r(f-z/4- ^7 



'r(z/i)r(z/2)r(z/3)r(z/4)r {n-vi-v^- v^) v {12-^2-^4- z/7) ' 



(D.2) 
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K^in; 0, z/2, z/g, 0, z/g, 1^6, = i 7r"(77' 



2— 2n n / 2\n— Si/j 



r(Ei/j - n) 



r(i/2)r(i/3)r(i/6)r (n - 1/2 - - ^^7) r (n - i/3 - - ^^e) r (f - Ei/,) 



,(D.3) 



2— 2n_n /'„2\n— Si-'i 



n 



/S:3(^; 0, z/2, 0, Z/4, Z/5, Z/6, Z/7) = i'-V{p 

^(f "^s) r(t~^6) r(^Z/5 + Z/6-f ) r(n-Z/2-Z/5-I^6-^'7) T (n - Z/4 - 1/5 - z/e - ^^7) 

r(z/2)r(z/4)r(z/5)r(z/6)r (n - 1/5 - ^^e) r (f - Ez/,) ' 

In the following formula we assume that z/7 = — s is a non-positive integer: 



(D.4) 



\2—2n„nl^2\n—YjV, 



TT 



r(Ei/i - n) 



r(f +5-1/2) r(f-i/3) r(f +S-Z/6) r(n+s-i/i-i/2-i/6) 
r(i/2)r(i/3)r(i/6)r (f - ei/,) r(n+s- 1/2-1/5) 

r (z/2 + z/6 - S - f ) -S, f - Z/3, z/2 + z/6 



3-^2 



5-2 



(D.5) 



r (z/i + z/2 + z/6 - S - 

Here, 3F2 denotes a generalized hypergeometric scries. In fact, we have a terminating 3F2 
series of unit argument, since one of the upper parameters is equal to —s. This may be 
considered just as a compact representation of a finite sum containing (s + 1) terms, each 
term being a product of F functions. 

Using the results for the boundary integrals with z/7 > 0, one can also calculate the 
planar forward-scattering double-box diagram (sec Fig. 3b), 



K,{n- 1, 1, 1, 1, 1, 1, 1) = -7r^-^^(-p2)-^- V 



■ + ■ 



16 TT^ 



-44-7r'+24C3+0(£) 



(D.6) 



and similar integrals with higher powers of the propagators. 



D.2: K2 integrals 

In the following two formulae, we also assume that 1/7 



-s is a non-positive integer: 



K2{n] 0, z/2, z/3, z/4, z/5, z/6, -s) = i 



2 — 2n „n / „2 \ n— Si/, 



TT 



F(Sz/i - n) 



r(f-z/2-^'4) r(f-z/3-z/5) F^f-z/e) F(n-z/2-z/3-z/4-z/6)F(n+s-z/3-z/4-z/5-z/6) 

X 7 \ 



F(z/2)F(z/5)F(z/6)F - Ez/j) F (n- z/2 -1/4-1/5) F (n- 1/3 -1/5-1/6) 

F(i/3 + l/5 + l/6-f) 



-S, i/3, ^ - i/2 - z/4 



^F (i/3 + 1/5 + l/g - S - f ) ^ V n-l/2-i/4-i/6, i^3 + i^5 + i^6-S-f 



(D.7) 



K^in; Z/i, 0, z/3, z/4, z/5, z/6, -s) = i^-^-^-^p2^n-Eu, r (j.^^ _ ^) 

f(|-z/4 + s) f(|-z/6 + s) F(n-z/i -Z/3-Z/4-Z/6 +s) F(n-z/3-z/4 -Z/5-Z/6 +s) 

X p 

F(z/i)F(z/4)F(z/5)F(z/6)F (n- 1/4-1/6-^5) F (^-1/1-1/3-1/4-1/5-1/6 +s 

XF(z/4 + Z/6-.-|) 3F2f -^,1-^5,^4 + ^/6-.-! 
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(D.8) 



Here, we also have terminating 3F2 series of unit argument, containing (s + 1) terms. 
However, in some cases they may contain less terms. For example, for an integer > 
the number of terms in 3F2 from ()D.8|) is min(s + 1,2/5). In particular, when s = or 
z/5 = 1 we get just one term. 

When = 0, eq. ()D.7|1 also contains just one term, since one of the upper parameters 
in 3F2 vanishes, so that 3F2 = 1. Moreover, the corresponding result can be extended to 
an arbitrary value of z/7: 

K2{n; 0, z/2, 0, u,, u,, z/g, z/7) = i'^-^V (p^)-^-^ T (Ez/, - n) 

r(f -z/5) r(f -Z/g) r(f -Z/2-Z/4) r(z/5 + Z/6-f ) r(n-Z/4-Z/5-Z/6-Z/7) 

r(z/2)r(z/5)r(z/6)r(z/5+z/6+z/7-t) r (n - z/5 - z/e) r (f - Sz/,) 

An important special case of eq. ()D.7|1 is z/7 = (s = 0). In this case, the 3F2 function 
is equal to 1 (since one of the upper parameters is zero) and we get0 

K2{n; 0, Z/2, Z/3, Z/4, z/5, Z/g, 0) = i2-2n^n^^2)„^S.,: ^ (^^^ _ ^) 

r(f-Z/2-I^4) r(|-Z/3-Z/5) r^f-Z/e) r(r2-Z/2-Z/3-Z/4-Z/6) r(n - Z/3 - Z/4 - Z/5 - Z/g) 
X p r . 

r(i^2)r(z/5)r(z/6)r {n-U2-Ui-U6) r (n-us-u^-ue) T - Sz/jJ 

(D.IO) 

This result corresponds to the last diagram shown in Fig. 5. 



Appendix E: Two-loop results for the ghost-gluon 

vertex 

Here we present the unrenormalized expressions for two-loop contributions to the scalar 
functions occurring in the ghost-gluon vertex ()2.9|) . in all on-shell limits of interest. To 
calculate these functions, the same algorithms (and the same REDUCE program) as for 
the three-gluon vertex have been employed. The two-loop diagrams contributing to the 
ghost-gluon vertex are shown in Fig. 2 of ref . PH] • Their renormalization is similar to that 
of the three-gluon vertex (see section 8 and Appendix B of ref. fS]); the renormalization 
factor Zi should be used. 

E.l: Non-zero gluon momentum squared 

^^This diagram has been considered in ref. '40^, using the method of negative-dimensional integration. 
We note that in their result (11) (or in the definition (5)) the parameters m and n (corresponding, up 
to a sign, to 1^4 and 1^3 in our eq. (jD.10|) ') should be interchanged. Also, {~tt)^ should read vr^ (they 
use Euclidean metric). We are grateful to the authors of |4()j for confirming these misprints. The other 
diagrams considered in jlHI (see also 41 ) correspond to integrals that can be obtained by repeated use 
of one-loop formulae. The results are given in eqs. (18), (21) and (23) of jlO], and they correspond to 
the special cases of our eqs. (jP.SII (at s = 0), ljD.2p (at i/7 = 0) and IID.4p (at 1/7 = 0), respectively. 
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£V 864 288 16^ 1728^ 64^'^^ 72^ 384 ^ 
117 1 

~32 

57245 _ 895 2 _ 179 1 4 _ 9139 3 2t_81^t 
^ 5184 ~ 864^ ~ Is"^^ ~ 64^ ~ 10368^ ^ 32^ ^ ~ 32^^'^ 

^1 _|_ _5^f2 3—77-2^2 f2 _ _|_ 19 2c3 

1280 ^ 1728^ 384 ^ 16^^^ 64^ 768 ^ 

4 



e^\72 2A 18^ 192 ^ 576^ 384 ^ 32^ 
384 ^ 32^ 

1 /779 31 2 1^ 1093^ 1 2^ 1^^ 331^2 1 2^2 
+ sl432 + T44" +4^^ + W^+96^^+2^^^+ 864^ "64^^ 

- + Y^TT ^ + — Ca^ - 

24959 293 2 25^ 1 4 5371 11 2^ 35^^ 3 4^ 
+ ^ + 432" + 24^^ + 80^ + 5184^ " 96" ^ + 16^^^+ 320^ ^ 

5549 2 7 2r^2 '^A/-2 ^ 4/-2 "^"^^^3 ^ 2 /-3 ^ /■ /-3 

+ 2592^ +192^^ -4^^^ -320^^ +48^^ +16^^^ 



+ 



,4 ^2 



1 / 25 7 ^ 7 . 45 . 

H — \ £ C C 

£3 V 48 32^ 384^ 256^ 

+1 _ 1^2 ,5 ^^2. _ 1.2 _ ^^2.2 _ ^.3 

£2 \^ 72 24 4^ 32 ^ 72^ 128 ^ 256^ 

I 3 2t3 1 c-4 

+ 256^ ^ - 16^ 

1 / 2407 3 2 1. 409 1 2. 9 1331 2 

^£ V I32" ~ 8^ ~ 4^^ ^ "96"^ ~ 8^ ^ ~ 32^^^ ~ 1728^ 

9 2 j-2 ^ 2 77 7 2 /"3 ^ /-3 /"4 

+ 128^^ "32^-^^ "128^' + 256^^ +64^^^ " 4^ 
26449 29 2 1^ 1 4 4651 11 2. 74^ 

-^-16^ -4^^-80" +W^-48"^-16^^^- 640"^ 
4471 2 73 2.2 87^ ^2 3 4^2 17.3 79 2.3 3^ . 



5184^ 192 32^"^ 320 ^ 8 768 8 
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,4 ^2 ^1/3 1 1 

V48 ^ 16^ 192^ 128^ ) 
1 /119 3 2 17^ 1 2e 263 2 7 2^2 29 3 
^?V72"~32^ ^ 24^ ~ 32^ ^ ^ 1152^ ~ 384^ ^ ~ 64^ 

I 1 2(^:3 

192 ^ 32^ 

1 /248 1 2 3^ 31 1 2^ 27^^ 307^2 1 2^2 
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"16^^^ " 128^ + 192'' ^ + 32^^^ " 16^ 

24947 103 2 17> 1 4 523^ 11 2. 25^^ 13 4^ 
+"648""72"^ "T^^"32^ + 432^+96^^"¥^^^" 640^^ 
I 7553 2 _ 2.2 _ H/^ _ _J_^ic2 _ , _Lvr2f3 + 1^^^ ^3 

2592^ 384 ^ 4 640 ^ 32^ 128 ^ 64^^^ 

- —A + 0(s), (E.4) 
640 ^ 32^ J ^ ^' ^ ^ 

A ^2 



^ ' ' ^ ^ (47r)"^ ^ l^i V 32 64^ 64^ 128^ 



V 16 32^ 96^ 128^ J 



£3 



+1 _ 1^2 ,93 ^ 2, _ ^.2 _ 1 ..2 _ ^.3 

£2 V 12 48 32^ 16 ^ 144^ 32 ^ 128^ 

I 3 2(^3 lf4 

+ 128^ ^ - 8^ 

1 / 35 67 2 17, 1681^ 17 o. 27,, 1439,0 

+ - TT^ Cs H C ttV Cs^ r 

£ V288 72 16^ 192^ 96 ^ 32^^^ 864 ^ 

+ 64^^ -32^^^ "64^ +128^^ +32^^^ "2^ 
33431 _ 1795 2 _ 125 . _ 3 4 19117 5 2-^ _ 9 > ^ 

^ 1728 ~ l32"^ ~ '2A^^ ~ 64^ ^ 1152 ^ ~ 16^ ^ ~ 4^^^ 

17 40^,55 2,2 171^,2 3 4 2 17,3 , 79 2,3 

"640^^" 648^ +96^^ " ^^^^ "128^^ " + 384^^ 

1 / 235 5 2 11 , 4 ,2 
+ £r216 + 36^ +108^+9^ 



4963 77 2 11. 157, ,, 61,2! , , 
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(47r)' 

1 / 23 2 2 

+ s IT08 " 9^ 
313 67 o 



6 3 6 



103, 103,2 



27 



54 



648 108 



7 _ 827 827 2 
3^^ ~ "sT^ ^ 162^ 



+ 



1 

18 



0{e), 



11 



11 



C + — C 



(E.7) 



1/239 1 2 16 94 
+^ll08 + 6^ -y^+27^ 
^ 2387 11 o 116 . _ 72 



648 



pV2''^)(p2,o,o) 



-C.T 



4 2 

(47r)" 



-2e 



6 



7 

^2 V9 ^ 6^ 



17 
18 
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+- \ £ f 

£ V27 18^ 54^ 



162 18 



54 
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1 / 73 5 2 32^ 188 .2 

+ir72 + T8^ 

4801 44 ^ 5 . 232. 1456 



432 

E.2: Non-zero in-ghost momentum squared 



44 2 5, 232, ^, , 1456,21 ^ 



^e^\72 192^ 64^ ^ 384^ ^ ^ 128^ 384^ ^ 

1 /355 11 2 1^ 305, 11 2^ 5^, 187,2 
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768 ^ 64 

37027 _ 791 2_175^_ 7 4 _ 263 23 2^ 33 1 4 

^ 2592 ~ 864^ ~ Is"^^ ~ 320^ ~ 16^ ^ 96^ ^ ^ 16^^^ ^ 160^ ^ 
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1/1 73 ^ 11 



£3 V 6 192'' 128' 

,1 ("l^. A^2_323 1 2 1 2 1 2 2 , ^.3,1 2,3 

eHl44 48 144^ 32 ^ 128^ 96 ^ 256^ 384 ^ 
1 /13391 31 2 5 14383 1 2^ 75 2 
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384 



^2^2 25 2 1013 ^^2^3 _ 1^^3\ 
^ 64^^ 128^ 96 ^ 32^^^ / 



315815 17 2 47^ 1 4 376861 31 2, 35^^ 7 4^ 
+ ^18r - - y^^- 32^ - l036^^+ 48" ^ + 16^^^+ 64o" ^ 
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£2 V 72 32 576 
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E.3: Non-zero out-ghost momentum squared 
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Figure 1: Contributions to the gg —>■ gg and qq gg processes which involve the three- 
gluon vertex 



P3 /"3 «3 




Pi ai P2 



Figure 2: Notation used for gluon momenta {pi + P2 + P3 = 0), Lorentz indices, /ij, and 
colour indices, 
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(a) 



(b) 



Figure 3: The planar two- loop three-point diagram (a) and an auxiliary four-point func- 
tion (b) 
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